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ABSTRACT 

The control of a single input-single output plant with parameter 
uncertainty has been approached using a minimax technique. A linear 
time invariant controller results which requires only partial state 
feedback and is optimal in the sense that it minimizes a quadratic 
criterion involving tracking error, control, and parameter uncertainty. 
Bounded input-bounded output stability is guaranteed provided the transfer 
function has only left half plane zeros. If uncertainties are bounded, 
it is always possible to stabilize the system when sufficient control 
amplitude is available. The number of states required to generate the 
control is equal to the system order less the number of zeros. 
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I. INTRODUCTION 

Control of a dynamic single input-single output system having parameter 
uncertainty and unknown process nonlinearity has been considered using a 
minimax technique. Specifically, we intend to show that a controller can be 
designed using a nominal model to insure a satisfactory performance of the 
system inspite of ignorance of system parameters. The problem is posed with 
the additional constraint that the controller is linear and that it requires 
only partial state feedback. 

II. 2 , 3] 

Lyapunov type synthesis techniques may be applied to control this 

class of plants. The resulting controller is highly nonlinear. The limitation 

of this approach is not only the complexity of the controller structure, but 

also the lack of any insight as how to determine the control signal amplitude. 

It is sometimes possible to obtain a design for this class of plant using 

14 ] 

Stochastic Control Theory. The general problem here is to determine the 
least favourable distribution for the uncertain parameter vector £. Another 
alternate approach is to use minimax control rule^ . 

In this paper, the problem has been approached by minimizing and maximizing 
a quadratic performance index involving the tracking error, the control, and 
the "uncertainty signal". The resulting controller is linear and is optimal 
in the sense that it minimizes the above performance criteria. The number of 
states required to generate the control is equal to the system order less the 
number of zeros. Bounded input-bounded output stability is guaranteed provided 
the transfer function is of minimum-phase type. If the uncertainties are 
bounded, it is always possible to stabilize the system, irrespective of the 
location of open loop poles of the system, if sufficient control is available. 
These results also apply for systems with rather general nonlinearities that 
do not involve the control. It has also been shown that the tracking error 
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vector admits an upper bound and that the bound can be made arbitrarily 
small. 


II. MOTIVATION OF THE PROBLEM 

We consider the single-input-single output system described as in 
(l)-(3): 


X, (s) b s m + b s m V . . . . + b,8 + b» .. 

1 = m m-1 1 0 N(s) 

U(s) n n-1 , , D(s) 

s + a .s + ....+ a. s + a_ 
n-1 1 0 


-m<n 


( 1 ) 


where X^(s) and U(s) represents the Laplace transforms of output X-^(t) and 
input u(t) respectively. The equation (1) may be reproduced in the state 
variable form^ 
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Since x^ + ^=x^, i = l»2, n-m, it can be seen from (1) that 


N(s) 
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X(s ) _ -1 1_ 

U(s) " (sI ~ A) h. D (s) 


n-m- 
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N (s) 
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N n (s) 


( 2 ) 
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where N (s), k=0,l,,..,m, are polynomials in s. Consider now the use of 

U” K 

a linear feedback law 

T 

u = - k. [k.x, + k„x~ + .... + k x ] = - k„k x ( 6 ) 

0 11 2 2 n-m n-m (F 

T 

where k = [k. ,k- , k , 0, ....0] is a constant vector, 

— l z n-m 

The eigenvalues of the closed loop system are then solutions of 

0 = | si - A + k () hk T | = | sI-A J | I + k Q (sI-A) -1 ]* 1 ! 

= | sI-A| (1 + k^ 1 (sl-A) _1 h) (5) 

where | A | is the determinate of matrix A. 

The last equality of (5) is obtained using the identity 

| 1 + C D j = 1 + D T C (6) 

£, D being vectors of compatible dimensions. Furthermore, combining (4) and 

the definition of k with (5) yields 

! 

T N(s)(k 1 + k.s + . . . . + k s n ~ m ) 1 

|sI-A + k o^ 1= °( s > I 1 + k 0 1 Hd(J) 1 

(7) 

- D < 8 > l 1 + k 0 n D(s? (e) 1 - ° 

(7) can also be obtained using Figure 1 which illustrates the system with 
feedback. The characteristic equation of the closed loop system is 

1 - Loop Gain = 1 + k^ 1 (sI-A)~ 1 h = 0, 

which agrees with (5) and (7). We know as kg-*=°, that zeros of (7) approach 

the n-1 finite zeros of N(s)»k(s) and one zero at -<*>. Hence if N(s), k(s) are 

Hurwitz polynomials with b , k >0, the system is stable for k n sufficiently 

I m n-m u 

large - regardless of the zeroB of D(s) = |sI-A|. Furthermore, response to any 
bounded input R(s) will be bounded. The problem is to choose the nonzero 
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elements of k so that in addition to stability the system exhibits hehavior 
which is in some sense good. Furthermore we should like a design procedure 
which can yield this good behavior despite uncertainties in the system parameters. 


III. SYSTEM DESCRIPTION, DEVELOPMENT AND PROBLEM FORMULATION 
Define a stable model 


+ $ rt r = 
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( 8 ) 


and a system 


x (n) + a x (n_1) + . 
1 n-1 1 


. . . + a^x, = b u m + b ,u 
0 1 m m-1 


m-1 


+ . . .+ b^u + bgU + 


f(x 1 ,...,x n _ m ,t) (9) 

where and r are the model output and reference input respectively, and 

f(x,, x„. . .x , t) is a nonlinear function, 
l z n— m 

If 

(1) b >p>0 and b s m + b . s m * + . . . +b,s + b n Hurwitz, and if 

m— m m-1 1 0 

(2) f(x^,X 2 ..., x nwm **) is a bounded continuous nonlinear function, 
we wish to show that the error 

e = (x , x, . . ., x ) T - 2 . “ x “ X (10) 

— 12 n-m j - — j - 

can be bounded with an arbitrarily small bound, despite imperfect knowledge of 

a i , b k , i = 0, 1 k = 0, 1, . . . m. This will be achieved by a 

linear feedback law which requires only partial state feedback and is in a 
sense optimal. Furthermore u will be similar in form to (4) with k(s) Hurwitz. 
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Remark 1. 


a) n,m, are integers which may be unknown but the difference (n-m) is 
assumed to be known. Also a lower bound of the highest order derivative 
coefficient of u is available. 

b) Note that (9) is identical to the linear system (1) except for 
addition of the nonlinear term, 

c) When a regulator is being designed it is admissable to let r=2=0. 

We proceed by rewriting (9) 


(n) + a , x, ( "- 1) 
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(j) 
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- I 

j=0 


a.x 


(j) 


j 1 


+ f(x. ... ., x t) 
l n— m 


(ID 


Integrating each side m times, gives 

(n-m) + a (n-m 1) + < < . + a x = g u + £(t) (12) 

1 n-m-1 1 01 


where £(t) is the m fold integral of the right side of (11) together with 
initial condition terms. In state variable form ,(12) may be written 
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500. 


= AqX + J3u +£(t) , 

* 

where x corresponds to the first n-m elements of x as in (10) 

, * r * (n-m-1) ,T , -.T 

i.e. x = l x 2 * • •> Xj^ '] = [Xj^, x 2 ,..., x n _ m J . 


(13) 


( 14 ) 
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Consequently subtracting ( 8 ) fjrpja (13) yields 

• 3 

1 = (x - £) = A Q e + ji (u + £) , u = (u- ~ r) ,£ = j £. (15) 

If 

■ * 

£(t), of course, generates differences between the model and x*/ £ ° 0 = u, e(t) 

approaches zero asymptotically since Aq has negative eigenvalues. Furthermore, 

even if u = 0, £ 0, e^ will be bounded for £(t) bounded. The problem is to 

a 

realize a u which will not only retain stability but allow/bound on _e to be 
made arbitrarily small. 


IV. CONTROLLER DESIGN 


For the purpose of controller design consider the criterion 


1 [ t '2 ~2 

J = j (e Qe + u R - £ L) dt 


(16) 


with R, L > 0 and Q positive definite, u, £ will be chosen to minimize and 

17] * * 

maximize J respectively. It is readily shown that the optimum u, £ are given 


by 


* IT 'IT 
u = - i 1 Pe, £ = f l A Pe 


(17) 


R — " L 

where P is the symmetric matrix satisfying 

PA q + aJp + Q - PS£ T P (| - - 0 . (18) 

1 T 

Furthermore, (i) Aq - — 66 P has negative eigenvalues 

(ii) P is the unique positive definite steady state 

solution of -P « PA q + aJp + Q - P 66 T P (^ - 7 ) , R<L. 

* * 

(iii) u, £ ] = min max J[eQ; u,£ ] = max min Jl^; u» £ ] 

u £ £ u 

1 T 

= ± e p e 

2 ^0 -^0 


(17) will, in one sense, yield a conservative design since £ is assumed to act 
in the most perverse manner. In another sense, however, £ need not abide 
by the rules of the game and may be using a smaller L than "agreed upon". This 
may require that R in turn be decreased. 



Remark 2. 


It should be noted that a unique positive definite solution P of (18) 

[81 

exists 1 J if R<L. In case, R<L the implication is that we are trying to balance 
the effect of the "uncertainty signal" with a larger amplitude of control 
signal than what is needed with R = L. Consequently for much of the rest of 
the paper, we shall restrict ourselves to the case R = L. In case R<L, 
it is still possible to show that proposition 1 (mentioned below) is true. 

In case R>L, nonuniqueness of solution P of (18) poses some problems in the 
subsequence analysis of the error bound and will be discussed in a future 
paper. 

For convenience, define as the last row or column of P. Then 

P-8.^ = l-R^ » £2 * ' * 1^,0,. • *»8] = ®£ n _ m - % • 

from (17), u is given by 

u ■ ~ |r e ■ - |k T e . (19*) 

Accordingly 

(j(s) = -fk T U,s,. . ., s 11 '” 1 " 1 ) 1 E(s)= - |k(s)E(s) (1%) 

where 


k(s) = 2 ^^ [l,s,s 2 . 


n-m-l,T 
3 ] 


Now it will be helpful, at this stage, to establish an important property 
of this control signal u that is outlined in the following proposition: 
Proposition 1 


With Aq stable and as defined in (13), k(s) is a stable polynomial for 

R=L. 


Proof : 

Consider the system defined by (15) and represented in Figure 2. The 
open loop transfer function, z(s), between £(s) and u(s) t is given by 

'<■) ■ I !? * e 2 ^ 


( 20 ) 
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where D^(s) = |sI-Aq| is a stable polynomial and k(s) is as defined in (19b). 
Now using (18) , we obtain 

2Re[z(s)]=z(s) + z(s) = ~ [^(sI-aJ)" 1 P.B + B T P (sI-Aq)' 1 ^ 

= |i T (sI-Aq) _1 [P(s+s) + P6 T (R -1 -L _1 ) BP +Q] (sI-Aq)" 1 ^ 

= 2 {Res} [jP (sI-Aq) _1 P(sI-A 0 )" 1 1] + B T (sI-aJ) - 1 Q(sl-A) _1 3 

for R = L. 

•k 

Here s denotes the complex conjugate of s. Since P, Q are positive definite, 
Re[z(s)J is nonnegative • for Re(s)>0. Therefore the transfer 

[9] 

function z(s) is positive real which implies from (20) 
that k(s) is Hurwitz. 

V. DERIVATION OF A BOUND ON SHE ERROR 
In order to determine a bound on je when u satisfies (17) but £(t) is 
arbitrary, let 

V(e) - ie T Pe (21) 

where P is the positive definite matrix satisfying (18). £ must now satisfy 

the differential equation (15) with u as in (19a); that is 

e = [A q -| B T p] e + 6 C(t). (22) 

The time derivative of (21) is 

V(e) = fe T lA Q - | M Tp ) T P + P (A q - B T P)e] e + £ (t)B T Pe 

= - |e T Qe - \ 6 2 (k T e) 2 (^ + £) + 5 (t)B (k T e) (23) 

where (18),(JL9a), and (22) have been used to refine the result. Clearly if 
|5(t)| is bounded ,(23) will be negative for | 1^1 | sufficiently large and admit 
an upper bound ||e||. First we must examine €(t), however, by considering 
the behavior of the full system described by (9), (10), (21). We shall return 
to complete our examination of (23) . 
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VI. STABILITY OF THE OVERALL SYSTEM 

Turning now to the total system, (9) may be written in state form, 
basically as in (2), 

x = Ax + hu + f (x, t) 


( 25 ) 


with X (x»t) = 10, 0 , . . ., 0 , f(x,t)] T . 

If we now define 


* 

rXi 


x = [-], A = 
it 


y 


(26) 


Z = [ y] , V - (▼ 1> v 2 . 


V 

m 


and combined (25), (8) and (10), to eliminate x, the result is 

Z “ I | r + f (e+j, t) 


Z = AZ + hu + 


A * “ A 0 




•* AZ + hu + v 


(27) 


where 


n— m n-m , 

v ■ 10,0,. ... 0, £ a y - B„r, 0, 0, ., f(£ + e, t) - l ,y,] 

i-0 1=1 


(n-m) entry 


(28) 


is bounded if the model input r, is bounded as assumed. If the 
"optimal" u given by (19a)is used in (27), the closed loop system satisfies 

B, 


Z = [A-|hk ] Z + h-r^r + v 

K — D ~ 


= IA - | h k T ] Z + v ’ 


(29) 


with 


k a [k. , k„ , ... k » 0,0, . . . 0] 

— x z n-m 

B n . 

v'= v + h — r = [0, 0, ... 0, v* ... v']. 
— — — p n— m n 


( 30 ) 
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The characteristic equation of the closed loop system is obtained as 
in (5)-(7), i.e. 

0= |s!-A+|hk X | = |sI-A| |l + (sI-A)" 1 | h k T | = |sI-A| (1 + | k T (sI-A) _1 h) 


- D(s) 11 + 1 bw <k i + V + • • • • + 'w (n ' m ' 1)l 


- D(s) + | N(s) k(s) (31) 

Now as R-H), (n-1) roots of (31) approach zeros of N(s)k(s), and the last root 
goes to infinity along the negative real axis. Since N(s) , k(s) are "stable 
polynomials", then (30) is stable for R adequately small. 

Consequently as the penalty R on control is reduced, permitting larger 
control amplitudes, the system (1) - (2) or (27) is stable for the feedback 
law (19). provided N(S) is a stable polynomial with b m *»0. 

Now since (29) yields a bounded Z with a bounded input v'. , all the elements 
of Z are bounded. In fact, bounds on v/ do not depend on Z but are determined 
mainly on From (29) - (31) it can be seen, in a straightforward manner, that 

+ 1 n . 

£(e.)=E,(s)= a I f^(s) v'(s) i=l,2, . . . , (n-m-1) 

1 1 D(s) + N(s)k(s) j=n-m 3 3 


£(v. )=V. (s/) g 

1 * D(s) + |N(s) fc(s) 

where f^(s), g^s) and (s) 
j j J 

order < (n-1). 


I Ig^Cs) + i- ^ (s) j v*(s), i=l,2, . , . ,m, 

j=n-m 3 3 3 

are polynomials in s, independent of R, and of 


As R-+0, |E i (s) |-> 0, i=l,2 ,..., (n-m-1) 


and 

|V i (s)| 


I ^(8) v'(s) 

1=n-m 3 3 

N(s)k(s) 


i=l,2,..,m. Thus an ultimate 
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bound on exists and is independent of R as R-K). Furthermore the error 
bound, i.e. the bound on | |ej | goes to zero as R-K). We explore this further, 

VII. FURTHER RESULTS ON THE BOUND OF THE ERROR 

Let us now examine 5(t). From (22) 

n-m oZ T 

5 (t) = 65 (t) = +J^ a i _ 1 e i + r k e 


n-m 


= e 


n-m 


+ , l (a i-l + R k i )e i 


i=l 


From (29) , on the otherhand 

n-m 

Z = e =v 1 -~ b k e + 7 a, ,y. - B^r + b — i 

n-m n-m 1 R m _ l-li 0 mB 


Thus 


n-m 


5(t) -V + Z + Vl’l 1 + 6 o ( f _1) r + t (s ' b » >i ‘ i 

1=1 


= ?1 + | (B-b m )k T e (32) 

where 

h ' v i + J 1 u i-i e i + t, i-i y i I + Ve 2 _1) 

has an upper bound which is independent of R and exists as R-»0. 

Returning to (23) with (32) replacing 5(t) 

V (e) - - i e T Qe - (kV^ ' H + fl> + h^ 

< " \ /Qe - \ (k T e) 2 B 2 (| + 

1 - 1 ^ - \ < k V s 2 i + r> + i«iLJ A I < 33 > 


- 2 


1 T- .1 
7 e Qe + - 


^ 1 ^ max 

<r4> 62 


W(e) 
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The last inequality above is obtained by maximizing the last two previous 

X 

terms with respect to | k _e | . It is possible to show in a straight forward 
manner that W(e) and hence V(e) is negative for all 


1= <e T e) 1/2 


Si I 


R = 
a 


1 'max 


6 


V 


^Vmin^L + 


where (^q^ m £ n denotes the minimum eigenvalue of Q. This, of course, is con- 
servative since we have permitted e to both maximize certain terms in (33) 
and to maximize | |e| | subject to W(e) >_ 0. 

The bound on | |e| | can now be found in standard fashion; 

(1) Determine V = maximum V(e) subject to | |ej | R 

U 3 

(2) Determine R^ = maximum j |e[ | subject to V(e) <_ 


e must asymptotically approach the region defined by V(e) <_ V^. 

R^ represents an ultimate bound on | |e| | since e^ asymptotically approaches 
the region defined by V(e) _< Vq but V(e) may be indefinite therein. Thus 
| | e | | is bounded by R^ where 



*/ ( Vma’x 
* ( Vmin 


- I? 


l'max 


*/ ^'Vmax 

* ^P^min^Q^min (^ + 7)6 

K L 


as shown in the Figure 4, This bound can be made arbitrarily small by relaxing 

the penalty R on the control u since l^jj max a dmits on upper bound as can be 

seen from previous section and (32). 

VII EXAMPLE 


To illustrate the preceeding analysis, consider the following open- 

loop unstable plant described by 

2 

x + x - x + x - Sinx + x = u+2u + u 


( 34 ) 
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wit h an arbitrary first order model described by 


y + 2y = r 

where r is a step input. It is to be noted that the minimum order of the 
model is specified by the difference between the number of poles and zeros 
of the plant. This is, in this example, one. With little manipulation 
as shown in the analysis and integrating (34) twice, the error equation can 
be expressed as 

e + 2e = £ + u (35) 

u = u-r 

* 

It is desired to find an optimal control u by maximizing w.r.t. £ and 
minimizing w.r.t. u the performance criterion 


J = 


(4e 2 + Ru 2 - L 5 2 ) dt, R=L 


subject to (35) , , 


(36) 


The resulting control is given by 

$ = - R _1 e (37) 

and is applied to the original system (34). 

The outputs of the plant and model are shown in Figure (5 ) for different 
values of R. It can be seen that the error decreases monotonically with 
the decrease of penalty R. Also plotted are the remaining states 
: x ^ ~ x, Xg= x for different values of R, . These are bounded as 

can be seen in figure (8) and (9). The control signal characteristics are 
shown in Figure (7), 
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IX. CONCLUSION 

A linear time invariant controller has been designed for a single input- 
single output system with parameter uncertainty. The number of states 
required to generate the control signal is equal to the system order less the 
number of zeros. This is obtained by minimizing a quadratic performance 
index involving the tracking error, the control signal and the "uncertainty 
signal". This however, yields a conservative design since the "uncertainty 
signal" is assumed to act in the most unfavourable manner. Bounded input- 
bounded output stability is guaranteed provided the transfer function is 
of minimum-phase type. If the uncertainties are bounded, it has been shown 
that the system can always be stabilized if sufficient control amplitude is 
available. These results also apply for systems with rather general nonlinear- 
ities that do not involve the control. We have also shown that the tracking 
error admits an upper bound and that the bound can be made arbitrarily small. 
Extension of this result to multivariable systems is not immediate and is 
currently under investigation. 
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